function z = func(x, y)
    z = 100.0 * (y - x.^2).^2 + (1 - x).^2;
end